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We present the next-to-leading order perturbative QCD prediction to the four-jet angular distri- 
butions used by experimental collaborations at LEP for measuring the QCD color charge factors. 
We compare our results to ALEPH data corrected to parton level. We perform a leading order 
"measurement" of the QCD color factor ratios by fitting the leading order perturbative predictions 
to the next-to- leading order result. Our result shows that in an experimental analysis for measuring 
the color charge factors the use of the Oiai,) QCD predictions instead of the 0(0^) results may shift 
the center of the fit by a relative factor of 1 + 2as in the Tr/Cf direction. 



I. INTRODUCTION 



In the first phase of operation of the Large Electron Positron Collider (LEP) four-jet events were primarily used 
for measuring the eigenvalues of the Casimir operators of the underlying symmetry group, the QCD color factors 
[0-p|. The values of these color charges test whether the dynamics is indeed described by an SU(3) symmetry. The 
\jr\ _ dependence of jet cross sections on the adjoint color charge appear at O(ck^). Several test variables with perturbative 
OO ' expansion starting at ©(a^) — so called four-jet angular distributions — were proposed as candidate observables 
C^ ■ with particular sensitivity to the gauge structure of the theory |g,g-g]. For a long time, however, perturbative 
CN QCD prediction for these variables at 0{a'l) had not been available, therefore the absolute normalization of the 

perturbative prediction could not be fixed. In order to circumvent this problem the experimental collaborations either 
fitted the strong coupling as well, or used normalized angular distributions of four jet events that were expected to be 
insensitive to renormalization scale dependence. The small scale dependence however, is but an indication and not a 
i-G ' proof of negligible radiative corrections: in principle the shape of the distribution can change from ©(a^) to O(a^). 
^~-i. Therefore, it is desirable to check explicitly the effect of the next-to-leading order corrections on these normalized 
O i' angle distributions. 

D , Recently the next-to-leading order corrections to various four-jet observables have been calculated |p|-pd|. These 

i-Q works depend crucially on the matrix elements for the relevant QCD subprocesses, i.e. for the e"'"e~ — > qqgg and 

^ , e^e^ — > qqQQ processes at one loop and for the e+e^ -^ qqggg and e+e~ — > qqQQg processes at tree level. The loop 

results became available due to the effort of two groups. In refs. [Ijjl^ Campbell, Glover and Miller made FORTRAN 

programs for the next-to-leading order squared matrix elements of the e^e~ — > 7* — + qqQQ and qqgg processes publicly 

S . available. In refs. ||lj,|l5[ Bern, Dixon, Kosower and Wienzierl gave analytic formulae for the helicity amplitudes of 

" " ' the same processes with the e'^e^ —> Z'^ —> four partons channel included as well. The helicity amplitudes for the 

five-parton processes have been known for a long time ||l^. Using the helicity amplitudes in refs. [|l4|-p^, Dixon and 

Signer calculated the next-to- leading order corrections to four-jet fractions for various clustering algorithms |S[|, as 

well as to the xbz angle distributio n |lO| . In previous publications we calculated several four-jet shape variables at 

the next-to- leading order accuracy [|1 1[| . In this paper we calculate the radiative corrections to the distributions of 

the commonly used angular shape variables ^ksw, ^nri '^34j ^^'^ repeat the calculation for the distribution of xbz- 

We use the matrix elements of refs. |lj,|l3l for the loop corrections, while calculated the helicity amplitudes of the 

relevant tree-level processes ourselves. 

Knowing the next-to-leading order corrections to these angular distributions, one would like to quantify their effect 
on the measurement of the QCD color charges. We estimate the systematic error coming from the use of leading order 
results instead of the next-to-leading order one in the fits for the color charge ratios x = Ca/Cf and y — Tn/Cp in 
the following way. We assume that next-to-leading order QCD is the true theory that describes the data. We fit the 
leading order prediction of the angular distributions with x and y left free to our next-to-leading order QCD results 
and determine these charge ratios from this fit. The central value of the "measured" charge ratios differs from the 
SU(3) values x = 9/4 and y = 3/8. This shift is the systematic bias that comes from the use of leading order fits in 
experimental analyses instead of the next-to-leading order ones. 



The rest of the paper is organized as follows. In section II we outline the structure of the numerical calculation and 
describe how we parametrise our results. In section III we present the complete 0(a|) predictions for the four standard 
angular distributions using two different jet algorithms: the Durham algorithm [|7| and the Cambridge algorithm 
proposed recently ||l^ . In section IV we perform the leading order fit of the color charges to our next-to- leading order 
results. Section V contains our conclusions. 

II. THE STRUCTURE OF THE NUMERICAL CALCULATION 

It is well known that the next-to-leading order correction is a sum of two integrals — the real and virtual corrections 
— that are separately divergent (in the infrared) in d = 4 dimensions. For infrared safe observables, for instance 
the four-jet angular distributions used in this work, their sum is finite. In order to obtain this finite correction, we 
use a slightly modified version of the dipole method of Catani and Seymour ||l^ that exposes the cancellation of the 
infrared singularities directly at the integrand level. The formal result of this cancellation is that the next-to-leading 
order correction is a sum of two finite integrals, 

NLO _ / j„NLO , / j„NLO 



a^^^^ = / da^^^^ + / dar"" , (1) 

where the first term is an integral over the available five-parton phase space (as defined by the jet observable) and 
the second one is an integral over the available four-parton phase space. 

Once the phase space integrations in eq. (|l|) are carried out, the next-to- leading order differential cross section for 
the four-jet observable O4 takes the general form 



(To dOi \ 2n J \ 27r 



2 /„ /.A^ \ 3 r o ..2 



1 da_^^^^^/aM^ ^^^(^^, , fa.{t^)CF 



BoA04)-^lTi— + CoA04 



(2) 



In this equation (Jq denotes the Born cross section for the process e+e^ -^ qq, s is the total cm. energy squared, /i is 
the renormalization scale, while -B04 and C04 are scale independent functions, Bd is the Born approximation and 
C04 is the radiative correction. We use the two-loop expression for the running coupling, 

with 

f . , ^ OsiMz). (Mz\ ,.. 

w{fi) = i-Po — ^ — l~rj ' ^' 

Po = Y^A - ItrNi , (5) 

01 ^Y^ A- 2CfTrNi - y^aTrNi , (6) 

with the normalization Tr = 1/2 in Tr(rT^'') — TuS""^. The numerical values presented in this letter were obtained 
at the Z° peak with Mz = 91.I87GeV, Tz = 2.49GeV, sin^ 6I = 0.23, a^iMz) = 0.118 and TVf = 5 light quark 
flavors. 

The Born approximation and the higher order correction are linear and quadratic forms of ratios of the color charges 

B4 = Bo + B^x + Byy, (7) 

and 

C4 = Co + CxX + Cyy + Cz z + Cxx x^ + C^y xy + Cyy y"^ . (8) 

At next-to-leading order the ratio z appears that is related to the square of a cubic Casimir, 



C3 = ^ Tr(r°T''rt'=)Tr(T^'=T''r'') , (9) 

a,fe,c— 1 



C3 



• ic ^ — ^ p3-. The Born functions Bi are obtained by integrating the fuUy exclusive 0{ag) ERT matrix elements |2l|] 

and were used by the experimental collaborations Ma. In the next section we present the Ci correction functions 
for the four different angular distributions. 

III. RESULTS 

In order to define the angular variables we denote the three-momenta of the four jets by pi, [i = 1, 2, 3, 4) and label 
jets in order of descending jet energy, such that jet 1 has the highest energy and jet 4 has the smallest. The four 
variables are 

1. the Korner-Schicrholz-Willrodt variable [|J,coS(/)ksw is the cosine of the average of two angles between planes 
spanned by the jets, 



0KSW = ;j 



(pi X Pi) ■ {p2 X P3)\ ( {pi X ps) ■ {p2 X Pi) 

arccos —— ttt^ ti- + arccos ' 



\pi X pa\\p2 X PsI / V |Pi ^ P?>\\'P2 X Pa 



(10) 



2. the modified Nachtmann-Reiter variable M, | cos6'^j^| is the absolute value of the cosine of the angle between 
the vectors pi — p2 and p^ — P4 , 

r.^.a* (Pi -P2) ■ (P3 -Pi) .,,^ 

\pi -P2\\pz-Pi\ 

3. cos 0^34 PI, the cosine of the angle between the two smallest energy jets, 

P3-P4 ,_,„. 

cos q;34 = I ^ 1 1 ^ I ; (12) 

b3|b4| 

4. the Bengtsson-Zerwas correlation |8J, | cosxbzI is the absolute value of the cosine of the angle between the plane 
spanned by jets 1 and 2 and that by jets 3 and 4, 

(pi XP2) • (P3 xp4) .-,„^ 

cosxBz = -n — ^TT^ — ^^ ; (13) 

blXp2|b3Xp4| 

We tabulate the numerical value of the next-to-leading order kinematic functions for these angular variables in the 
Appendix in Tables III- VI for the Durham clustering algorithm and in Tables VII-X for the Cambridge algorithm 
that was proposed recently |18| . Using this new algorithm, the hadronization corrections are expected to be much 
smaller therefore, the perturbative prediction is more reliable. The values in the tables were obtained by selecting 
four-jet events at a fixed jet resolution parameter ^cut = 0.008 which is the value used by the ALEPH collaboration 
[^. We do not show the value of the C^ functions because they turn out to be negligible. The C4 values were obtained 
according to eq. (||) with SU(3) values for the color charge ratios, x = 9/4, y = 3/8. Comparing the size of the 
corrections for these two algorithms, we see that in general the C,; functions in the case of the Cambridge algorithm 
are 10-20 % smaller. 

We use the numerical values for the kinematic functions to calculate the next-to-leading order QCD predictions 
for the SU(3) values x = 9/4, y = 3/8 according to eq. (pi) at x^ = pij \fs = 1. We compare our predictions for the 
Durham algorithm (solid histograms) to ALEPH data (diamonds) in Figs. 1-4. In order to make this comparison we 
normalize the histograms to one, therefore 

F(A ^ i5j(.) (14) 

in the plots. The qualitative agreement between data and theory is very good. Also shown in these figures our results 
for the Cambridge algorithm (dotted histograms). The statistical error of the Monte Carlo integrals is below 1.5% 



for the Durham algorithm and below 2 % in the case of the Cambridge algorithm in each of the bins, 
figures, the windows show polinoniial fits to the K factors of the normalized distributions defined as 



In the same 



K{z) 



1 dcTNLO 

(Tnlo d2; 



(^) 



1 dgLO 
7L0 dz 



(z) 



(15) 



where ctnlo = " + o is the next-to-leading order cross section. The x^/^doi of these fits is between (1.5-7)/20. 
The K factors for the | cos xbz| distributions, for the cos 034 distribution with Durham algorithm and for the | cos 0^^ \ 
distribution with the Cambridge algorithm are approximately constant 1 over the whole range, therefore the shape of 
the leading and next-to-leading order distributions are very similar in these cases. 
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FIG. 1. Comparison of the next-to- leading order QCD prediction for the cos^ksw distribution obtained using Durham 
(solid) and Cambridge (dotted) jet algorithm with ALEPH data (diamonds) . In the window the K factor of the distribution 
with Durham (solid) and with the Cambridge (dotted) algorithm is shown. 
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FIG. 2. Comparison of the next-to-leading order QCD prediction for the | cos 0f^^\ distribution obtained using Durham (solid) 
and Cambridge (dotted) jet algorithm with ALEPH data (diamonds). In the window the K factor of the distribution with 
Durham (solid) and with the Cambridge (dotted) algorithm is shown. 
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FIG. 3. Comparison of the next-to-leading order QCD prediction for the cos Q34 distribution obtained using Durham (sohd) 
and Cambridge (dotted) jet algorithm with ALEPH data (diamonds). In the window the K factor of the distribution with 
Durham (solid) and with the Cambridge (dotted) algorithm is shown. 
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FIG. 4. Comparison of the next-to-leading order QCD prediction for the |cosxbz| distribution obtained using Durham 
(solid) and Cambridge (dotted) jet algorithm with ALEPH data (diamonds). In the window the K factor of the distribution 
with Durham (solid) and with the Cambridge (dotted) algorithm is shown. 



IV. LEADING ORDER VERSUS NEXT-TO-LEADING ORDER 



A quantitative comparison of the data for the angular distributions to the next-to-leading order prediction de- 
composed in a quadratic form of the color factor ratios with group independent kincmatical functions as coefficients 
makes possible a simultaneous fit of the strong coupling and the color charge ratios. That procedure would require 
a full experimental analysis which is not our goal in the present paper. What we would like to achieve is to give 
a reliable estimate of the systematic theoretical uncertainty coming from the use of the leading order perturbative 
prediction instead of the next-to-leading order one in a color charge measurement. To this end we pretend that the 
ncxt-to-lcading order perturbative QCD is the "true" theory that describes data perfectly. We "produce" data using 
our next-to-leading order prediction with SU(3) values x = 9/4 and y — 3/8 and perform a leading order fit of x and 
y using the Bqj Bx and By functions. We use y^ minimalization to obtain the best values with 



2_\p 1 I Bo{zi) + xBx{zi) + yByjZj) 1 dcTNLO . A ^„. 

^ w} V CTo + xa^ + j/o-j, CTNLO dz ' ) 

where Wi is the statistical error of the norraaUzcd next-to-lcading order distribution in the ith bin, the summation 
runs over the bins and Oj(zi) (j = 0, a;, y, or NLO) defined as 

CTj = dzBj{z), (TNLO = / dz ^'"" (z) (17) 

As a check of the fit we also performed a linear fit to the not normalized distributions in the form 

X^^Y.^('^(^o{z^)+xB4z^) + yBy{z,))-^^^{z,)) , (18) 

where Wi is the statistical error of the next-to-leading order distribution in the ith bin, and with rj = [asCp / {'^t^))'^ 
fitted as well. The two procedures give the same result for x and y to very good accuracy. 

We performed the fit for each angular distribution separately, as well as for the four angular variables combined. 
Table I contains the results of these fits. We see that the shifts in the x-y values are quite large. Looking at the 
errors, one finds that the shift is significant only if the K factor of the corresponding distribution (see Figs. 1-4) 
is not constant 1. For those cases when the shapes of the leading order and the next-to- leading order distributions 
are very similar, i.e. the A" ~ 1, then the fits give values compatible with the canonical QCD values. The origin 
of the large errors in some fits is the global correlation between the two parameters x and y in the fit. In these 
cases — |cos0Jjpj|, cos 034 and |cosxbz| distributions — one cannot fit both variables reliably. Instead, one can 
fit either the ratio of the two parameters, or fix one parameter to the SU(3) value and fix the other. For instance, 
fixing X = 9/4 one obtains the fitted values for y as given in Table II. We observe from Figs. 1-4 and Table II 
that in those cases, when iiT ~ 1 the result of the fit is in agreement with SU(3) — Durham cos 034, |cosxbz| 
and Cambridge |cos0J,j^| distributions — , while for the rest of the distributions we obtain fit parameter differ- 
ent from the SU(3) value because the shapes of the leading and next-to-leading order distributions are different. 



TABLE I. 
correlations. 


Leading 


order fit of the color charge 


ratios to the next-to-leading 


order differential distributions 


of the 


angular 


Observable 






X 






y 


Durham algorithm 



cosc^KSw 2.21 ±0.05 0.58 ±0.07 

|cos(9nr| 1.41 ±1.43 0.08 ±0.11 

cos a34 2.08 ± 0.21 0.57 ± 0.23 

|cosxBz| 1.15 ±1.43 0.12 ±0.31 

all four 2.32 ± 0.03 0.29 ± 0.02 

Cambridge algorithm 

cos <;/>KSw 2.30 ± 0.08 0.52 ± 0.09 

Icos^nrI 0.99 ±2.70 0.21 ±0.31 

cos a34 0.34 ± 0.42 2.65 ± 0.48 

I cos XBz I 3.53 ± 2.80 0.82 ± 0.68 

all four 2.29 ± 0.05 0.45 ± 0.03 



TABLE II. Leading order fit the color charge ratio y to the next-to-leading order differential distributions of the angular 
correlations with x = 9/4 fixed. 

Observable Durham algorithm Cambridge algorithm 

cos <;/)KSw 0.57 ±0.06 0.55 ± 0.08 

Icos^nrI 0.15 ±0.03 0.36 ±0.05 

cos 034 0.39 ± 0.05 0.56 ± 0.08 

I cos XBZ I 0.35 ± 0.05 0.51 ± 0.06 

all four 0.31 ± 0.02 0.46 ± 0.03 



We show the result of the combined fit of all four variables in Fig. 5 in the form of 68.3 % and 95 % confidence level 
contours in the x-y plane with ellipses centered on the best x-y pair. There are five contours sitting on three different 



centers in each plot. The fits with both 1- and 2-cr contours were obtained using all four angular distributions with 
all bins included. The fit with only l-cr contour shown corresponds to the "ALEPH choice" : using all four variables 
with fit ranges 0.1 < | cosxbzI, | cos0^j^| < 0.9 and —0.8 < cosa34, cos^ksw < 0.8. 

We observe from Fig. 5., that the leading order fit results in overestimating the Ca/Cf ratio by 2-3% no matter 
which clustering algorithm is used. For the TbJCf ratio the leading order fit underestimates the result by 20-30% 
of a next-to-leading order fit when the the Durham algorithm is used, while in the case of Cambridge clustering 
the leading order fit gives an overestimate of about 20%. This systematic bias appears significant in both cases. 
Although the two parameters are slightly correlated when all four variables are used, the fit is reliable. The result 
of the fit depends on the jet algorithm because the different jet finders lead to different jet momenta from which 
our test variables are built. We also see that constraining the fit range as the ALEPH collaboration did does not 
alter our conclusions significantly. We would like to emphasize that the significant shift from the SU(3) values does 
not mean the exclusion of QCD, but simply gives an estimate of the systematic theoretical error in the color charge 
measurements when leading order fits are used. 
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FIG. 5. Confidence level contours of the leading order fits of the color charges x and y for the Durham and Cambridge 
clustering at ycut ~ 0.008. 



One may ask how the light gluino exclusion significance changes in the recent analysis of Csikor and Fodor [p2J, 
which used the results of four-jet analyses, if one takes into account the systematic theoretical error discussed above. 
Assuming that the shifts of x and y are similar in the light gluino extension of QCD, our conclusion suggests that 
the radiative corrections induce a shift of order 2as times the tree-level value for x and y. Lacking this piece of 



information Csikor and Fodor have increased the axes of the error elhpses by a factor of as times the theoretical x 
and y values. Implementing our results to a Csikor-Fodor type analysis for the four-jet events would decrease their 
confidence levels for the light gluino exclusion from 99.9% (Csikor-Fodor value) to ~98%, which is, however, still 
much higher than a 2-cr exclusion. 

V. CONCLUSION 

In this paper we presented the next-to-leading order corrections to the group independent kinematical functions 
of the four standard four-jet angular distributions, cos^ksWj I cos0Jj^|, cosa34, and | cos^bzI with jets defined with 
two — the Durham and the Cambridge — clustering algorithms at j/cut = 0.008. These results were obtained 
with a general purpose Monte Carlo program called DEBRECEN [g3| that can be used to calculate the differential 
distribution of any other four-jet quantity at the next-to-leading order accuracy in electron-positron annihilation. The 
I cos xbz I distribution using the Durham clustering algorithm was first calculated by Signer [ [L0[ . Our result agrees 
with his within statistical errors. The results for the other three distributions with Durham clustering as well as for 
all distributions when the Cambridge algorithm is used are new. We compared our results to data obtained by the 
ALEPH collaboration corrected to parton level and found very good qualitative agreement. We have also presented 
the K factors of the distributions using both jet clustering algorithms. 

Having the ncxt-to-leading order perturbative QCD prediction at our disposal we made an estimate of the systematic 
theoretical error of the QCD color charge measurements due to the use of leading order group independent kinematical 
functions. We found that the use of the 0(as) QCD predictions instead of the 0(0;^) results may shift the center of 
the fit by a relative factor of about 1 -I- 2as in the Tj^/Cp direction, while the best Ca/Cf value is hardly affected. 
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Scientific Research Fund grant OTKA T-016613 and the Research Group in Physics of the Hungarian Academy of 
Sciences, Debrecen. 
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VI. APPENDIX 



TABLE 

used. 


III. Next-to-leading 


order kinematical functions to the 


cos 4>KBw angular 


distribution. 


The Durham jet algorithm is 


cos (?!)KSW 


Ci 


Co 


a 


Cy 


a. 


Cxy 


Cyy 


-0.950 


1190.6 ±25.0 


-159.3 ±6.7 


655.2 ± 8.9 - 


-1504.8 ±11.8 


129.7 ±3.6 


-235.6 ±5.2 


-128.0 ±0.7 


-0.850 


686.5 ±21.2 


-86.7 ±5.8 


381.8 ± 7.6 


-889.8 ±9.3 


70.7 ±2.9 


-109.6 ±3.8 


-127.9 ±1.0 


-0.750 


696.9 ± 19.4 


-67.5 ±5.1 


353.1 ± 7.2 


-826.6 ±9.8 


73.3 ± 2.8 


-82.4 ±3.6 


-157.4 ±1.4 


-0.650 


707.5 ± 22.6 


-73.8 ±5.3 


352.2 ± 7.8 


-822.8 ±9.7 


73.9 ± 3.0 


-61.5 ±3.8 


-179.4 ±1.7 


-0.550 


722.0 ± 20.9 


-73.7 ±5.0 


366.5 ± 8.1 


-851.9 ±9.9 


70.9 ± 2.7 


-48.2 ±3.9 


-196.7 ±1.8 


-0.450 


716.2 ±19.7 


-77.9 ± 5.2 


373.9 ± 8.2 


-881.7 ±11.0 


67.0 ± 2.4 


-32.6 ±3.9 


-203.8 ±2.1 


-0.350 


708.5 ± 19.6 


-87.8 ±5.9 


390.4 ± 7.9 


-893.5 ±10.4 


57.8 ±2.3 


-10.3 ±3.5 


-206.5 ±2.0 


-0.250 


763.1 ±20.6 


-78.2 ±6.5 


418.1 ±9.2 


-959.0 ±11.7 


57.5 ±2.4 


-2.0 ±3.3 


-207.1 ±1.9 


-0.150 


752.5 ±19.7 


-77.7 ±14.0 


435.4 ± 10.6 


-982.8 ±11.7 


46.9 ±2.2 


10.4 ±3.0 


-197.2 ±1.8 


-0.050 


730.7 ± 18.8 


-104.0 ±7.3 


457.3 ± 9.7 - 


-1039.0 ±12.7 


41.1 ±2.1 


16.1 ±2.6 


-191.6 ±2.0 


0.050 


665.5 ±18.7 


-90.7 ±7.3 


420.7 ± 9.9 


-934.2 ±12.6 


31.6 ±2.1 


30.0 ±2.5 


-176.1 ±1.8 


0.150 


652.9 ±18.0 


-84.7 ±7.2 


431.0 ±9.6 


-935.4 ±13.3 


21.2 ±2.1 


40.6 ± 2.6 


-157.4 ±1.5 


0.250 


631.6 ±17.0 


-86.2 ±7.5 


424.2 ± 9.8 


-904.9 ± 13.4 


17.5 ±2.1 


40.2 ±2.1 


-147.9 ±1.4 


0.350 


633.7 ±15.6 


-83.3 ±6.7 


433.1 ±9.5 


-929.5 ±13.4 


14.5 ±2.2 


44.2 ±2.2 


-137.6 ±1.3 


0.450 


576.9 ± 29.9 


-105.1 ±7.3 


404.6 ± 42.6 


-923.4 ±12.7 


18.9 ±13.3 


47.5 ±2.2 


-127.2 ±1.2 


0.550 


628.8 ±29.5 


-100.3 ±8.5 


484.4 ± 42.6 


-980.6 ±13.5 


-4.0 ±13.5 


52.6 ±2.6 


-122.8 ±1.1 


0.650 


628.4 ±16.2 


-114.0 ±8.7 


483.6 ± 12.3 - 


-1041.5 ±14.2 


2.3 ±3.3 


59.3 ±2.9 


-119.7 ±1.0 


0.750 


711.1 ±15.9 


-132.8 ±9.8 


578.1 ± 13.9 - 


-1197.2 ±15.5 


-9.8 ±3.7 


68.2 ±3.1 


-118.2 ±0.9 


0.850 


836.8 ±15.6 - 


-169.5 ±10.8 


675.8 ± 14.2 - 


-1444.1 ± 17.4 


-3.3 ±4.1 


74.0 ±3.7 


-130.3 ±1.1 


0.950 


1820.0 ± 22.9 


-399.8 ±12.9 


1547.2 ± 16.8 - 


-3258.7 ±23.5 


-29.1 ±4.8 


168.9 ±5.1 


-252.6 ±1.3 



TABLE IV. Next-to- leading order kinematical functions to the | cos &nrI angular distribution. The Durham jet algorithm is 
used. 



jcose^Niii 


Ci 


Co 


C. 


Cy 


C 


Oxy 


'-'1/1/ 


0.025 


1184.8 ±26.6 


-164.6 ±10.5 


757.6 ± 16.8 


-1769.9 ±24.6 


46.2 ±5.5 


161.7 ±7.8 


-459.4 ± 4.3 


0.075 


1130.0 ±35.3 


-170.6 ±12.3 


757.8 ±21.3 


-1755.2 ±26.9 


36.7 ±6.7 


156.8 ±8.6 


-450.0 ± 4.4 


0.125 


1204.4 ± 37.0 


-164.8 ±12.0 


780.8 ± 22.2 


-1783.0 ±25.5 


43.0 ±7.1 


150.8 ±8.2 


-451.7 ±4.1 


0.175 


1190.0 ± 39.9 


-170.9 ±14.7 


764.6 ± 20.2 


-1813.6 ±28.1 


52.0 ±6.6 


142.6 ±8.5 


-451.4 ±4.1 


0.225 


1239.0 ± 37.9 


-172.1 ±13.9 


801.1 ±18.7 


-1828.6 ±28.7 


47.1 ±5.9 


139.5 ± 7.5 


-432.3 ± 4.0 


0.275 


1247.5 ± 32.2 


-172.4 ±12.4 


796.1 ± 18.3 


-1832.1 ±25.7 


54.2 ±5.8 


115.2 ±8.0 


-423.6 ± 4.0 


0.325 


1268.2 ±33.5 


-189.6 ±15.3 


828.9 ±19.7 


-1885.8 ±25.7 


54.7 ±6.4 


100.1 ±8.5 


-408.9 ±3.5 


0.375 


1311.9 ±40.1 


-187.4 ±14.6 


824.8 ± 20.1 


-1903.0 ±26.0 


68.5 ±6.3 


74.3 ± 7.7 


-386.8 ±3.2 


0.425 


1358.3 ±41.6 


-177.3 ±13.0 


888.3 ± 20.6 


-1986.6 ±26.8 


56.5 ±6.3 


55.9 ±7.4 


-374.3 ±3.2 


0.475 


1456.3 ± 38.0 


-201.0 ±13.5 


916.8 ±20.5 


-2044.7 ± 27.4 


74.9 ± 5.9 


38.5 ± 7.4 


-356.7 ±3.0 


0.525 


1443.4 ± 41.3 


-225.8 ±14.3 


951.0 ±20.6 


-2130.4 ± 26.3 


69.0 ± 5.8 


34.1 ±6.9 


-342.0 ±3.1 


0.575 


1570.0 ± 60.5 


-208.6 ± 14.3 


980.2 ± 29.9 


-2198.9 ±26.5 


89.4 ±6.3 


-11.4 ±7.2 


-317.0 ±2.5 


0.625 


1662.7 ±59.7 


-249.2 ± 14.6 


1085.3 ± 28.5 


-2305.5 ± 26.3 


79.1 ±5.8 


-30.2 ±6.7 


-300.8 ±2.5 


0.675 


1806.9 ± 60.9 


-233.2 ±15.9 


1146.4 ± 30.3 


-2448.1 ±31.4 


92.6 ±6.1 


-61.2 ±7.1 


-276.8 ±2.2 


0.725 


1689.3 ±62.1 


-279.2 ±18.5 


1129.5 ±30.7 


-2503.3 ± 30.8 


93.3 ± 6.3 


-83.3 ±6.8 


-249.5 ±1.8 


0.775 


1913.3 ±50.4 


-263.6 ±16.6 


1244.8 ± 23.8 


-2605.8 ± 29.0 


93.6 ± 6.9 


-105.4 ±6.2 


-226.0 ±1.6 


0.825 


1931.5 ±52.0 


-290.4 ±17.1 


1269.4 ± 22.1 


-2670.3 ± 28.0 


101.2 ±6.2 


-139.6 ±6.7 


-200.1 ±1.4 


0.875 


1907.6 ± 50.7 


-284.7 ±17.9 


1254.5 ±22.7 


-2691.0 ±30.5 


106.3 ±5.7 


-160.8 ±7.6 


-173.4 ±1.3 


0.925 


1979.0 ± 56.1 


-302.0 ±31.4 


1297.6 ±26.1 


-2777.3 ±29.1 


116.5 ±6.3 


-196.4 ±7.4 


-151.9 ±1.1 


0.975 


2426.4 ± 68.3 


-399.1 ±19.4 


1658.0 ± 27.8 


-3469.3 ± 33.7 


122.2 ±6.5 


-241.3 ±8.3 


-137.8 ±0.8 
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TABLE V. Next-to-leading 
used. 


order kinematical functions to the cos 034 angular 


distribution. 


The Durham jet algorithm is 


cos 034 


Ca 


Co 


a 


Cy 


^xx 


Cxy 


Cyy 


-0.950 


1038.4 ± 17.9 


-227.3 ±9.1 


901.2 ±12.1 


-1888.6 ±17.5 


-20.6 ±2.8 


77.1 ±2.9 


-106.0 ±0.8 


-0.850 


988.9 ± 19.5 


-187.7 ±9.0 


830.9 ± 12.8 


-1717.5 ± 17.9 


-19.0 ±3.0 


73.8 ±2.9 


-109.0 ±0.8 


-0.750 


948.7 ± 19.4 


-177.0 ±9.2 


777.1 ± 13.9 


-1628.8 ±16.6 


-10.4 ±3.5 


67.2 ±2.7 


-116.9 ±0.9 


-0.650 


911.1 ±22.2 


-162.4 ±9.6 


716.5 ±13.0 


-1500.5 ±15.4 


-0.9 ±2.9 


54.8 ±2.7 


-123.8 ±1.0 


-0.550 


899.6 ±26.1 


-144.1 ±7.9 


684.6 ± 14.7 


-1427.9 ±15.6 


3.4 ±3.1 


47.9 ± 2.7 


-132.5 ±1.0 


-0.450 


849.9 ± 22.1 


-131.0 ±15.0 


629.7 ±14.3 


-1381.1 ± 15.8 


12.9 ±3.4 


43.0 ±3.0 


-141.0 ±1.1 


-0.350 


827.5 ±21.1 


-133.5 ±7.3 


597.8 ±11.5 


-1268.9 ±15.7 


16.0 ±3.2 


38.1 ±2.9 


-148.8 ±1.2 


-0.250 


816.0 ±19.0 


-126.8 ±6.8 


566.3 ±11.6 


-1235.8 ±14.7 


25.9 ±2.9 


26.5 ±2.8 


-158.2 ±1.2 


-0.150 


810.4 ±19.8 


-125.3 ±7.0 


546.5 ± 12.2 


-1200.7 ±14.2 


32.2 ±3.6 


20.0 ±3.2 


-171.5 ±1.3 


-0.050 


843.6 ±21.4 


-118.9 ±7.4 


534.6 ±13.5 


-1187.5 ±14.9 


44.4 ±3.9 


7.4 ±3.2 


-187.6 ±1.5 


0.050 


829.9 ±21.7 


-126.0 ±6.7 


510.9 ±10.3 


-1152.5 ±13.3 


53.0 ±2.7 


-1.4 ±3.4 


-203.9 ± 1.6 


0.150 


879.4 ± 23.4 


-105.2 ±7.2 


484.1 ± 10.1 


-1100.1 ± 13.4 


69.4 ±2.9 


-14.5 ±3.9 


-219.1 ±1.8 


0.250 


834.2 ±21.6 


-104.4 ±6.2 


450.9 ± 8.9 


-1055.2 ± 12.4 


74.0 ± 2.7 


-28.7 ±4.2 


-233.1 ±2.1 


0.350 


837.2 ± 23.9 


-93.4 ±5.6 


426.1 ±8.7 


-999.6 ±11.9 


83.0 ±3.1 


-43.0 ±4.5 


-244.0 ± 2.0 


0.450 


828.0 ± 24.0 


-79.6 ± 5.3 


382.7 ±7.8 


-915.3 ±11.1 


93.2 ±3.4 


-55.2 ±4.9 


-249.9 ±2.1 


0.550 


793.7 ± 22.2 


-72.9 ±4.5 


343.9 ± 8.9 


-854.5 ± 10.3 


100.1 ±2.7 


-70.2 ±5.0 


-252.9 ±2.1 


0.650 


719.8 ±21.9 


-59.1 ±5.0 


305.2 ± 8.4 


-763.9 ± 8.4 


94.0 ± 2.9 


-74.7 ±5.0 


-238.4 ±2.1 


0.750 


546.0 ± 15.0 


-52.5 ±4.1 


238.4 ± 5.2 


-602.4 ± 7.5 


73.0 ±2.5 


-66.1 ±4.0 


-184.8 ±2.1 


0.850 


245.2 ± 9.3 


-24.0 ± 2.9 


130.1 ±3.4 


-302.8 ±5.2 


24.8 ±1.7 


-31.7 ±2.4 


-62.8 ±0.9 


0.950 


12.5 ±2.4 


-1.9 ±1.0 


9.2 ±1.3 


-17.7 ±1.4 


0.2 ±0.2 


-0.4 ±0.2 


-0.9 ±0.1 



TABLE VL Noxt-to- leading order kinematical functions to the | cosxbzI angular distribution. The Durham jet algorithm is 
used. 



JCOSXBZ 


C4 


Co 


Cx 


Cy 


Cxx 


<-yxy 


'-'1/1/ 


0.025 


1133.9 ±39.7 


-171.8 ± 14.1 


748.6 ± 24.7 


-1623.5 ±30.9 


29.2 ±6.3 


161.8 ±8.2 


-400.9 ± 4.4 


0.075 


1087.1 ±47.9 


-170.2 ±17.3 


724.3 ± 35.3 


-1689.3 ±31.6 


34.9 ± 9.9 


164.9 ±9.2 


-395.9 ±4.2 


0.125 


1147.7 ±49.5 


-163.0 ±18.2 


771.0 ±35.7 


-1640.7 ±26.6 


23.1 ±11.0 


153.2 ±7.8 


-395.1 ±4.1 


0.175 


1075.4 ± 45.9 


-160.2 ±16.6 


719.1 ±27.9 


-1665.6 ±26.8 


33.6 ±8.9 


159.0 ±9.3 


-391.4 ±4.2 


0.225 


1144.5 ±42.5 


-177.8 ±16.9 


752.8 ± 27.3 


-1659.9 ±25.7 


36.9 ±7.5 


135.4 ±8.9 


-381.4 ±3.9 


0.275 


1110.3 ±43.5 


-172.6 ±17.9 


741.9 ± 29.0 


-1634.2 ±27.2 


33.4 ±8.2 


121.7 ±8.1 


-361.1 ±3.7 


0.325 


1152.2 ±42.3 


-142.2 ± 18.6 


744.5 ± 28.2 


-1651.3 ±26.3 


39.0 ± 7.8 


110.2 ±7.7 


-350.3 ±3.1 


0.375 


1181.2 ±43.1 


-175.4 ±16.9 


794.6 ± 28.2 


-1741.2 ±26.3 


35.9 ± 7.3 


106.1 ±6.8 


-340.5 ±3.3 


0.425 


1210.8 ±47.8 


-155.7 ±16.7 


798.5 ± 36.3 


-1693.3 ±26.2 


37.9 ±10.1 


72.5 ± 7.5 


-330.6 ±3.2 


0.475 


1213.1 ±44.5 


-185.8 ±17.7 


762.7 ± 34.5 


-1749.9 ±25.3 


66.2 ±9.4 


56.7 ±7.7 


-322.3 ± 2.9 


0.525 


1236.6 ± 38.9 


-195.1 ±18.2 


819.6 ± 24.7 


-1801.7 ±23.2 


53.0 ±6.6 


44.5 ±6.3 


-310.7 ±2.8 


0.575 


1337.2 ± 42.9 


-182.4 ± 18.0 


847.8 ± 24.9 


-1855.7 ±24.3 


65.6 ±6.7 


19.1 ±5.8 


-299.8 ± 2.5 


0.625 


1383.5 ± 44.4 


-201.5 ±16.0 


891.7 ±23.7 


-1983.3 ±24.3 


72.4 ±6.5 


-2.5 ±5.8 


-292.6 ± 2.3 


0.675 


1415.8 ±44.2 


-197.9 ±15.9 


916.5 ±24.2 


-2050.8 ± 25.1 


74.8 ± 6.4 


-22.5 ±6.1 


-279.9 ±2.0 


0.725 


1551.9 ±44.0 


-283.2 ± 89.5 


1039.4 ± 77.4 


-2150.5 ±25.3 


74.9 ± 15.7 


-45.4 ±6.4 


-268.6 ±2.0 


0.775 


1659.4 ± 47.8 


-152.9 ±90.5 


1020.2 ± 77.7 


-2357.4 ± 26.3 


99.2 ±15.8 


-75.9 ±6.7 


-262.1 ±1.9 


0.825 


1828.2 ±48.2 


-251.4 ±17.7 


1166.8 ±26.1 


-2553.8 ± 25.9 


105.5 ±6.3 


-102.4 ±7.0 


-248.7 ±1.6 


0.875 


1791.5 ±245.1 


-316.9 ±32.0 


1236.4 ±91.5 


-2863.7 ±27.2 


110.1 ±8.8 


-146.7 ±7.0 


-243.1 ± 1.6 


0.925 


2455.5 ± 246.3 


-318.9 ±33.3 


1528.1 ±91.8 


-3230.3 ± 28.2 


151.2 ±9.1 


-217.0 ±7.9 


-249.3 ± 1.5 


0.975 


4804.5 ± 70.2 


-731.4 ±36.2 


3109.1 ±34.0 


-6806.4 ±41.7 


320.2 ±8.8 


-552.9 ± 12.4 


-445.9 ± 1.9 
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TABLE VII. Noxt-to-leading order kinematical functions to the cos <^ksw angu' 
is used. 


lar distribution. 


The Cambridge jet algorithm 


cos 0KSW 


Ca 


Co 


a 


Cy 


^xx 


^xy 


Cyy 


-0.950 


1085.9 ±44.0 


-287.0 ±8.1 


713.0 ±12.8 


-1754.9 ±17.4 


135.1 ±7.2 


-280.6 ±8.0 


-143.7 ±0.9 


-0.850 


658.4 ± 35.4 


-164.3 ±7.1 


404.4 ± 10.0 


-986.2 ±12.6 


81.7 ±4.9 


-134.9 ±6.1 


-134.0 ±1.2 


-0.750 


617.4 ±27.0 


-139.8 ±8.8 


366.6 ± 9.6 


-932.2 ±13.9 


78.5 ±3.7 


-108.1 ±6.0 


-170.1 ±1.6 


-0.650 


620.0 ± 23.2 


-124.2 ±8.8 


355.9 ±9.5 


-924.0 ±14.5 


74.7 ±3.1 


-71.1 ±5.7 


-196.1 ±2.0 


-0.550 


662.8 ±23.0 


-130.5 ±7.9 


377.7 ±10.5 


-961.9 ±13.0 


75.7 ±3.2 


-59.5 ±4.6 


-210.4 ±2.2 


-0.450 


605.2 ± 24.4 


-145.8 ±8.2 


380.8 ± 10.6 


-954.7 ±13.1 


60.8 ±3.0 


-30.6 ±4.4 


-216.4 ±2.5 


-0.350 


658.4 ±23.8 


-143.1 ±7.8 


405.4 ± 10.4 


-1004.1 ± 14.8 


60.9 ±2.7 


-13.3 ±4.4 


-224.1 ±2.7 


-0.250 


653.0 ±25.4 


-148.1 ±7.6 


421.2 ±11.6 


-1035.8 ±15.6 


54.3 ± 2.6 


-4.2 ±4.6 


-219.6 ±2.5 


-0.150 


610.7 ±25.0 


-173.7 ±8.8 


431.9 ±12.1 


-1044.3 ±16.1 


44.3 ± 2.8 


13.4 ±4.3 


-213.7 ±2.4 


-0.050 


603.3 ± 23.6 


-176.6 ±12.3 


449.2 ± 15.2 


-1078.7 ±16.5 


36.9 ± 3.4 


16.1 ±3.8 


-192.3 ±2.2 


0.050 


548.4 ±21.4 


-147.6 ±7.5 


401.7 ±10.6 


-970.0 ± 14.6 


30.8 ±2.5 


28.0 ±2.6 


-171.0 ±2.0 


0.150 


514.7 ±19.3 


-161.5 ±7.3 


414.3 ±10.2 


-952.0 ±13.2 


18.2 ±2.4 


37.2 ±2.9 


-161.0 ±1.8 


0.250 


508.4 ±17.1 


-153.8 ±6.9 


407.1 ± 10.0 


-942.3 ±12.7 


17.3 ±2.2 


38.5 ± 2.9 


-147.1 ±1.6 


0.350 


500.3 ± 17.0 


-158.0 ±7.4 


425.3 ± 13.3 


-974.9 ±14.6 


9.9 ±3.2 


42.6 ±2.2 


-140.6 ±1.5 


0.450 


497.7 ±18.5 


-164.4 ±7.5 


429.0 ± 13.9 


-996.6 ±13.9 


10.0 ±3.4 


44.7 ±2.4 


-129.7 ±1.3 


0.550 


471.1 ±16.9 


-177.1 ±7.9 


444.7 ±12.0 


-1027.1 ±15.3 


1.1 ±2.9 


52.3 ±2.7 


-124.8 ±1.3 


0.650 


510.0 ±16.8 


-195.1 ±7.8 


489.6 ±11.6 


-1085.2 ±15.0 


-4.2 ±3.1 


56.5 ±2.9 


-120.9 ±1.2 


0.750 


560.0 ± 17.4 


-207.1 ±8.3 


545.7 ±12.3 


-1217.8 ±15.2 


-8.4 ±3.2 


65.2 ±3.4 


-118.7 ±1.0 


0.850 


681.4 ±22.0 


-264.5 ±9.5 


687.5 ± 14.6 


-1514.9 ±18.8 


-16.0 ±3.7 


78.5 ±3.9 


-129.7 ±1.2 


0.950 


1493.1 ± 25.4 


-630.2 ±14.1 


1574.3 ± 18.8 


-3435.0 ± 24.6 


-48.4 ±5.7 


176.7 ±6.1 


-252.8 ± 1.4 


TABLE VIII. Next-to- leading order kinematical functions to the cos ^^r angular distribution, 
is used. 


The Cambridge jet algorithm 


Icose^NRJ 


Ca 


Co 


a 


Cy 


Cxx 


^xy 


^yy 


0.025 


1016.1 ±30.2 


-233.6 ± 14.4 


750.2 ± 20.8 


-1827.3 ±28.5 


34.8 ± 6.8 


160.4 ±9.5 


-449.3 ± 4.6 


0.075 


1023.7 ±38.2 


-237.2 ±12.6 


734.5 ± 22.5 


-1817.0 ±32.1 


46.1 ± 7.2 


136.2 ±10.2 


-450.5 ±4.6 


0.125 


1011.2 ±44.1 


-224.1 ± 13.7 


756.4 ± 23.2 


-1848.0 ±30.1 


31.7 ±7.3 


157.7 ±9.7 


-441.6 ±5.0 


0.175 


1038.5 ± 44.3 


-249.2 ± 13.3 


760.8 ± 21.9 


-1867.3 ±27.9 


44.8 ± 6.8 


126.8 ±10.1 


-450.3 ± 4.5 


0.225 


1037.3 ±42.3 


-240.7 ±11.0 


760.2 ± 24.2 


-1867.6 ±29.5 


43.9 ±6.7 


117.8 ±11.4 


-430.8 ± 4.0 


0.275 


1112.1 ±42.0 


-253.5 ±12.1 


807.4 ± 23.2 


-1937.9 ±27.2 


47.7 ± 6.8 


117.0 ± 11.3 


-426.1 ±4.6 


0.325 


1078.7 ±37.7 


-260.3 ±13.6 


808.9 ± 21.4 


-1937.7 ±27.0 


44.1 ±6.7 


99.2 ±9.0 


-416.5 ±4.0 


0.375 


1074.5 ± 37.2 


-277.2 ±13.6 


797.5 ±21.5 


-1997.0 ±29.2 


57.6 ±6.9 


87.5 ±10.1 


-403.2 ±3.8 


0.425 


1147.9 ±39.6 


-314.1 ± 15.2 


896.8 ± 22.6 


-2103.0 ±28.6 


46.6 ± 6.4 


57.1 ± 10.0 


-385.7 ±3.7 


0.475 


1222.8 ± 43.0 


-294.4 ± 13.3 


911.7 ±23.0 


-2166.9 ±30.0 


58.9 ±6.6 


38.1 ±7.9 


-370.6 ±3.9 


0.525 


1238.6 ± 46.5 


-343.9 ± 15.0 


944.9 ± 23.2 


-2233.5 ±27.5 


62.5 ±6.7 


29.5 ±8.0 


-352.4 ±3.7 


0.575 


1295.5 ±45.3 


-351.4 ± 15.1 


969.0 ± 23.0 


-2346.7 ± 30.8 


78.8 ± 6.0 


-4.6 ±7.1 


-334.3 ±3.5 


0.625 


1371.8 ±45.0 


-393.9 ± 14.4 


1069.7 ±22.8 


-2488.9 ±31.3 


72.3 ±6.0 


-39.0 ±8.7 


-315.0 ±3.3 


0.675 


1403.5 ± 55.5 


-452.3 ±17.6 


1134.1 ±25.8 


-2606.4 ± 34.4 


75.7 ±6.7 


-73.7 ±8.6 


-294.4 ±2.8 


0.725 


1504.7 ±61.1 


-455.9 ±25.1 


1183.5 ±33.8 


-2732.3 ±34.7 


85.6 ±8.9 


-86.4 ±8.1 


-274.2 ±2.8 


0.775 


1479.6 ± 59.4 


-539.4 ± 17.8 


1202.1 ±29.3 


-2802.7 ±33.0 


99.7 ±8.7 


-126.3 ±8.0 


-250.6 ± 2.5 


0.825 


1527.5 ±63.3 


-566.6 ±20.2 


1268.4 ±32.6 


-2896.6 ±34.1 


94.9 ±8.2 


-142.1 ±8.7 


-230.8 ± 2.4 


0.875 


1585.8 ± 64.0 


-629.1 ±21.1 


1318.0 ±27.1 


-2971.0 ±34.6 


110.0 ±7.3 


-196.0 ±10.0 


-204.7 ±2.0 


0.925 


1706.3 ±87.7 


-625.3 ±20.9 


1357.8 ± 29.4 


-3096.3 ± 37.4 


130.8 ± 10.4 


-234.3 ± 11.1 


-180.9 ±1.8 


0.975 


2243.9 ± 92.8 


-842.0 ± 23.9 


1818.8 ±32.5 


-4041.0 ±44.8 


160.3 ±13.9 


-329.7 ±15.2 


-171.8 ±1.5 
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TABLE IX. Next-to-leading order kinematical functions 
is used. 


to the cos 034 angular distribution. 


The Cambridge 


jet algorithm 


cos 034 


C4 


Co 


a 


Cy 


^XX 


^xy 


Cyy 


-0.950 


827.2 ±21.4 


-373.0 ±9.4 


904.6 ± 13.6 


-1956.5 ±16.9 


-30.4 ±3.3 


79.5 ±3.1 


-105.8 ±0.9 


-0.850 


786.6 ±41.6 


-333.1 ± 10.3 


819.7 ±20.6 


-1750.8 ±19.7 


-22.2 ±3.0 


70.2 ± 3.3 


-108.9 ±0.9 


-0.750 


731.5 ±40.2 


-310.3 ±11.6 


760.7 ± 19.6 


-1686.0 ±17.9 


-15.0 ±3.2 


64.5 ±3.2 


-118.2 ±1.0 


-0.650 


726.5 ±21.4 


-285.2 ±10.0 


718.2 ±13.6 


-1607.3 ±18.0 


-6.1 ±3.2 


56.4 ±3.1 


-126.8 ±1.1 


-0.550 


691.0 ±30.1 


-257.9 ±8.5 


666.9 ± 15.9 


-1473.9 ±18.3 


-4.3 ±4.2 


48.0 ± 3.3 


-133.0 ±1.2 


-0.450 


704.7 ± 32.9 


-242.1 ±9.5 


635.0 ±16.8 


-1439.9 ±18.6 


8.5 ±4.2 


40.2 ±2.8 


-143.0 ±1.2 


-0.350 


692.8 ± 23.8 


-211.5 ±9.5 


590.3 ±13.9 


-1344.2 ±21.1 


14.2 ±3.4 


35.6 ±5.2 


-149.0 ±1.3 


-0.250 


690.1 ±21.9 


-203.2 ± 7.9 


563.4 ± 12.8 


-1266.4 ±20.2 


20.0 ±3.1 


24.9 ± 5.2 


-158.9 ±1.4 


-0.150 


670.2 ± 22.9 


-193.3 ±8.8 


531.8 ±12.5 


-1250.5 ±17.0 


28.9 ±3.0 


17.4 ±3.4 


-174.4 ±1.5 


-0.050 


699.2 ± 23.9 


-180.2 ±8.3 


511.6 ±12.3 


-1212.6 ±16.9 


38.8 ±3.1 


13.9 ±3.7 


-185.9 ±1.9 


0.050 


720.4 ±30.1 


-176.5 ±7.3 


499.8 ± 13.8 


-1163.4 ±14.7 


48.0 ± 2.8 


-7.5 ±4.1 


-202.7 ±1.8 


0.150 


769.0 ±31.7 


-158.3 ±7.4 


487.1 ± 14.7 


-1185.0 ±16.3 


62.2 ±3.0 


-10.5 ±4.7 


-223.5 ±2.1 


0.250 


791.6 ±33.7 


-163.2 ±7.5 


477.7 ±11.2 


-1163.8 ±15.8 


74.3 ± 4.3 


-31.5 ±4.9 


-239.8 ±2.1 


0.350 


790.2 ± 28.0 


-159.0 ±6.8 


451.3 ±10.4 


-1119.0 ±14.8 


86.1 ±3.6 


-54.3 ± 5.6 


-260.3 ± 2.3 


0.450 


816.2 ±33.8 


-143.8 ±6.7 


426.7 ±11.2 


-1095.9 ±14.4 


100.2 ±4.2 


-69.3 ± 5.8 


-273.9 ± 2.9 


0.550 


733.4 ± 30.9 


-148.9 ±7.1 


370.3 ± 10.9 


-1033.3 ±14.7 


110.4 ±4.2 


-98.8 ± 7.0 


-277.3 ± 2.7 


0.650 


634.0 ± 35.4 


-143.2 ±6.3 


333.7 ±9.4 


-942.1 ± 14.4 


98.3 ± 5.9 


-97.3 ±7.1 


-263.6 ± 2.7 


0.750 


434.4 ± 22.9 


-122.8 ±10.4 


235.7 ±13.3 


-732.8 ±11.0 


80.7 ±4.8 


-91.5 ±5.9 


-201.8 ±2.4 


0.850 


150.1 ± 14.0 


-79.6 ± 3.6 


135.6 ±7.6 


-355.6 ±8.4 


20.2 ±3.0 


-41.5 ±4.7 


-69.4 ±1.4 


0.950 


0.9 ±1.7 


-7.2 ±1.0 


5.1 ±1.0 


-13.8 ±1.3 


0.5 ±0.2 


-0.9 ±0.2 


-0.7 ±0.1 



TABLE X. Noxt-to- leading order kinematical functions to the | cosxbzI angular distribution. The Cambridge jet algorithm 
is used. 



JCOSXBZ 


Ci 


Co 


Cx 


Cy 


Cxx 


Oxy 


^yy 


0.025 


950.6 ± 40.9 


-233.8 ± 16.6 


718.7 ±29.8 


-1751.2 ±31.7 


28.4 ±8.0 


156.9 ±8.1 


-398.6 ±4.5 


0.075 


968.5 ± 53.7 


-242.8 ± 16.5 


762.8 ±35.1 


-1695.1 ±30.3 


11.9 ±9.3 


151.7 ±8.6 


-395.7 ±4.9 


0.125 


949.3 ± 52.6 


-245.5 ±16.7 


718.9 ± 34.7 


-1784.5 ±34.0 


31.7 ±9.0 


168.0 ±9.6 


-404.9 ± 5.2 


0.175 


989.3 ± 45.2 


-246.1 ± 17.8 


751.7 ±32.2 


-1713.5 ±32.8 


22.8 ±9.2 


141.4 ±8.4 


-390.4 ±4.6 


0.225 


959.7 ±45.3 


-260.6 ±21.7 


722.5 ±35.1 


-1707.9 ±32.1 


34.8 ±9.7 


138.6 ±8.3 


-383.3 ±4.7 


0.275 


940.8 ± 45.5 


-288.5 ±38.1 


736.2 ± 38.3 


-1702.2 ±29.2 


33.1 ±8.6 


114.1 ±8.2 


-373.5 ±4.4 


0.325 


989.1 ±48.3 


-224.8 ± 39.8 


710.5 ±38.1 


-1744.9 ±29.7 


45.0 ±9.1 


108.5 ±8.7 


-367.0 ± 4.3 


0.375 


1011.7 ±52.2 


-257.5 ±17.5 


781.7 ±30.0 


-1804.6 ±31.9 


30.5 ±8.2 


96.4 ±9.3 


-358.6 ±4.0 


0.425 


1055.0 ±49.7 


-262.7 ±17.0 


767.8 ± 30.7 


-1870.0 ±31.5 


55.2 ±8.3 


75.1 ±9.0 


-345.1 ±3.3 


0.475 


1015.4 ±71.3 


-378.7 ±99.1 


862.4 ± 63.3 


-1902.1 ±32.4 


31.6 ±15.2 


63.9 ± 8.4 


-344.3 ± 3.8 


0.525 


1110.7 ±89.5 


-194.5 ±98.2 


764.8 ±68.1 


-1943.5 ±32.6 


63.6 ± 15.0 


46.9 ± 7.7 


-328.1 ±3.4 


0.575 


1136.4 ± 75.5 


-297.7 ±18.6 


834.4 ± 39.9 


-1991.1 ±29.0 


67.3 ± 8.0 


8.2 ±9.6 


-308.5 ±3.0 


0.625 


1156.8 ±50.0 


-313.0 ±18.6 


893.3 ± 27.0 


-2106.2 ±28.9 


58.5 ±7.7 


-7.5 ±9.7 


-311.5 ±3.5 


0.675 


1223.0 ± 60.9 


-379.1 ± 17.0 


947.1 ±33.5 


-2229.2 ±28.1 


74.0 ± 7.6 


-31.0 ±7.4 


-298.0 ±2.7 


0.725 


1257.6 ± 65.3 


-391.6 ±27.5 


1005.7 ± 40.8 


-2386.4 ± 30.4 


72.3 ±10.5 


-51.3 ±8.4 


-288.4 ±2.6 


0.775 


1355.8 ± 56.8 


-435.0 ± 27.6 


1058.5 ± 34.6 


-2530.1 ±28.3 


93.4 ±9.7 


-92.2 ±8.9 


-273.9 ± 2.4 


0.825 


1386.3 ± 70.3 


-527.7 ±21.2 


1179.4 ±29.8 


-2785.3 ±30.9 


87.8 ±8.0 


-119.5 ±8.4 


-269.5 ±2.2 


0.875 


1636.2 ± 70.6 


-568.2 ±21.3 


1311.0 ±31.0 


-2996.6 ±31.6 


109.3 ±8.7 


-165.5 ±9.1 


-256.2 ±2.0 


0.925 


1994.9 ±111.6 


-674.3 ±22.1 


1553.8 ± 63.4 


-3496.2 ± 38.8 


146.4 ±15.4 


-261.0 ± 11.2 


-264.5 ± 1.8 


0.975 


4032.8 ± 121.5 


-1362.2 ±25.5 


3169.4 ± 67.0 


-7444.4 ±49.1 


329.2 ± 18.8 


-646.6 ± 17.3 


-473.7 ±2.1 
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